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We analyze the thermodynamical properties, at finite density and nonzero temperature, of the 
(1 + l)-dimensionaI chiral Gross-Neveu model (the NJL2 model), using the exact inhomogeneous 
(crystalline) condensate solutions to the gap equation. The continuous chiral symmetry of the 
model plays a crucial role, and the thermodynamics leads to a broken phase with a periodic spiral 
condensate, the "chiral spiral" , as a thermodynamically preferred limit of the more general "twisted 
kink crystal" solution of the gap equation. This situation should be contrasted with the Gross-Neveu 
model, which has a discrete chiral symmetry, and for which the phase diagram has a crystalline phase 
with a periodic kink crystal. We use a combination of analytic, numerical and Ginzburg-Landau 
techniques to study various parts of the phase diagram. 



I. INTRODUCTION 



The phase diagram of interacting fermion systems at finite density and temperature is a general problem with appli- 
cations in a wide range of physical contexts. Well-studied examples include the Peierls-Frohlich model of conduction 
[l| , the Gorkov-Bogoliubov-de Gennes approach to superconductivity Q , and the Nambu- Jona Lasinio (NJL) model 
of symmetry breaking in particle physics l3|. Strongly interacting theories such as quantum chromodynamics (QCD) 
exhibit a rich phase diagram structure It is known that chiral symmetry plays a key role, and computationally 

the large Nf and large Nc limits must be addressed carefully 0, [I]- A (1 -I- l)-dimensional version of the NJL model, 
the NJL2 model [also known as the chiral Gross-Neveu model, XGN2] is of interest because it captures some important 
features of QCD, such as asymptotic freedom, dynamical mass generation, a large Nf limit, and the breaking of a 
continuous chiral symmetry IS , 
equation, found recently in [l3, 



10 , [11|, . In this paper we use the exact crystalline solutions to the associated gap 
14| , to study the temperature-density phase diagram of this NJL2 system. The result 



of our thermodynamical analysis confirms the physical picture proposed in [15| | that there is a phase transition at a 
critical temperature Tc from a massless phase to a broken phase with a helical condensate (the "chiral spiral"), of the 
complex Larkin-Ovchinikov-Fulde-Ferrell (LOFF) form. The resulting phase diagram [see below. Figure [S], is very 
different from that of the non-chiral Gross-Neveu (GN2) model, which has just a discrete, rather than continuous, 
chiral symmetry. In the GN2 model there is also a region of the phase diagram with a crystalline order parameter 
(l6t . but the structure of the phase diagram is very different [see below, Figure [7]. This crystalhne phase of GN2 



has been clearly seen in a recent lattice analysis jl7l |. extending an important earlier lattice analysis [18[. In this 



paper we explain in detail the role of the chiral symmetry [continuous versus discrete] in determining the form of 
the phase diagram. The chiral spiral phase of the NJL2 model has also been studied in the AdS/QCD framework 
PH . These one dimensional models are of course simplified models of more realistic (3 -|- l)-dimensional systems, 
but important lessons can still be learned concerning the appearance of crystalline structures in the phase diagram 
[igl. [20|. Furthermore, their solubility permits a detailed study of the relation between real and imaginary chemical 
potential [2^ . 

Our analysis is ultimately based on solving the gap equation for inhomogeneous condensates. Initially, the phase 
diagram of the N JL2 and GN2 models was studied assuming homogeneous condensates [2^ [l^l , but this assumption 
does not capture certain aspects of the true physical phase diagram [16,, 17, l^]. Of course, finding inhomogeneous 
solutions to the gap equation is a much more difficult technical problem, but the massless NJL2 and GN2 models have 
remarkable symmetry properties that enable one to find the general periodic condensate solutions [l^, [3] • This fact 
is due to a deep connection between the Bogoliubov-de Gennes effective hamiltonian of the NJL2 system, and certain 
one dimensional integrable hierarchies (25l. l26l. [27j . These exact solutions are characterized by a finite number of 
parameters, and to describe the phase diagram one must minimize the thermodynamical grand potential with respect 
to these parameters in order to determine the form of the condensate in a given region of the (T, ^) plane. This 
thermodynamical analysis is performed in this paper. 

In Section II we briefly review the analytical solution of the inhomogeneous gap equation. In Section III we identify 
the special role played by rescaling and phase rotation symmetries in the NJL2 model. The thermodynamics of the 
NJL2 model is discussed in terms of a spiral condensate in Section IV and in terms of the general twisted kink crystal 
in Section V. In Section VI we contrast this analysis with the case of the GN2 model, which has just a discrete chiral 
symmetry. In Section VII we apply a Ginzburg-Landau analysis to study the region of the phase diagrams of both 
the GN2 and NJL2 models, in the vicinity of the relevant "tricritical point" . We conclude with a summary of our 
results and a discussion of the implications for more complicated models. 
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II. SOLVING THE INHOMOGENEOUS GAP EQUATION 

The NJL2 model is described by the following (1 + l)-dimensional Lagrangian with both scalar and pseudoscalar 
four-fcrmion interaction terms: 



NJL 



(2.1) 



This model has a continuous chiral symmetry: -0 — > e*''' "■(/;. The GN2 model has just the scalar four-fermion 
interaction term: 
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(2.2) 



and has a discrete chiral symmetry: tp 7^-0. We study these models in the large Nf limit where the semiclassical 
approximation applies and chiral symmetry breaking can be studied [H, |2^ . 

By a Hubbard-Stratonovich transformation, the four-fermion interaction terms can be expressed in terms of scalar 
and pseudo-scalar bosonic condensate fields, S and If (respectively), which are conveniently expressed in terms of a 
complex condensate field: A E — ifl. For GN2 we only have S, and so the condensate field A is real. The general 
NJL2 system can be described equivalently by the effective Lagrangian: 



£ = V 



»^-^(l-7')A-i(l + 7')A* 



which is now quadratic in the fermion fields. The corresponding single particle fermionic Hamiltoninan is 



H 



dx 



1 



n'x:+7"U(l-7')A-7T(l + 7')A 



(2.3) 



(2.4) 



With the choice of the Dirac matrices as 70 = fi, 71 = —ia2 and 7^ = cr3, the Hamiltonian 

H 



takes the form: 



-^1 A(x) 



A*(x) 



"' dx 



(2.5) 



This Hamiltonian is also known as the Bogoliubov-de Gennes (BdG) (or Andreev) Hamiltonian in the superconduc- 
tivity literature 0, [Slj . 

There are two equivalent perspectives on studying the semiclassical gap equation for static condensates. The first, 
a Hartree-Fock approach, is to solve the single particle equation (the Bogoliubov-de Gennes equation) 



Hi) = Eij} 



(2.6) 



subject to the consistency condition relating the condensate field to the expectation values of the scalar and pseu- 
doscalar fermionic bilinears: 



(2.7) 



A second approach is to integrate out the fermionic field in (j2.3p and obtain an effective action (per fermion flavor) 
for the condensate field: 



1 



(fxlAl^ - ilndet 



i(l-7')A-^(l+7')A* 



The gap equation for the condensate field is obtained by looking for the stationary points of S'ofr[A] 

SScs 







6A* 



1 



2gmf 



A{x) 



5A{xy 



In det 



z^-i(l-7^)A(x)-i(l-f7')A(x)' 



(2.9) 



It is straightforward to solve this gap equation when the condensate field A is uniform, but it is more technically 
challenging to solve it for an inhomogeneous condensate field A(x). Nevertheless, in one spatial dimension it is possible 
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to find the most general bounded quasi-periodic solution to this gap equation [l^, [l^ . The general solution has the 
form of a "twisted kink crystal" , described below. 

A useful quantity for solving the inhomogeneous gap equation (|2.9[1 is the resolvent R{x;E), the coincident-point 
limit of the Gor'kov Green's function G{x,y;E) corresponding to the Hamiltonian (|2.5p 

R{x;E)^{x\j^\x) . (2.10) 

For a static condensate the gap equation (|2.9[) can be written as 

A(a;) = -iNfg^ tr d^e [7° (l + 7') R{x; E)] (2.11) 

The solution of the gap equation relies on the remarkable fact that in one spatial dimension the resolvent (itself a 
2x2 matrix) must satisfy a simple first order matrix differential equation 

l-Rix;E)a,=z[( J,, '^l^A^Ra,] (2.12) 



dx ' ' ' '~ \A*ix) -E 

This equation is known as the Eilenberger equation in the superconductivity literature [10, [3l| , and as the Dickey 
equation in mathematical physics [l^, [sJl • The Dickey-Eilenberger equation follows immediately from the fact that 
the resolvent can be written as a product of two linearly independent solutions: 

R{x-E) = ^{xl:i^l + ^2i^l)ai (2.13) 

where W = i{'4>i U2'4'2) is the Wronskian of two independent solutions ?/'i,2 of Hijj = Eip. 

The inhomogeneous gap equation (|2.9p can be solved by the following simple ansatz [l^ [l3| for the resolvent 

R(x-E^-M(E\( + h{E)A-zA'\ 

^^'''-^> ~^^^>[b{E)A* +iA*' a{E) + \A\-' ) ^^'^^^ 

where Af{E), a{E) and b{E) are functions of the energy E, and are to be determined. This particular ansatz is 
motivated by the gap equation (|2.1ip that relates the off-diagonal component of the resolvent with A. The ansatz 
(|2.14p automatically solves the diagonal part of the Eilenberger equation (|2.12p . while the off-diagonal part requires 
that the condensate field A satisfy the complex nonlinear Scrodinger equation (NLSE): 

A" - 2|A|2 A + i {b{E) - 2E) A' - 2 {a{E) ^ Eb{E)) A = (2.15) 

The advantage of this ansatz approach is that the NLSE (|2.15p can be solved in closed form, and its general solution 
has the form of a twisted kink crystal, described in detail in [l3| and summarized below in the next section. The 
associated energy functions JV{E), a{E) and b{E) are simple fimctions of E. For the NJL2 model, there is a further 
consistency condition required to satisfy the gap equation (|2.1ip : for an inhomogeneous condensate, the part of the 
off-diagonal resolvent proportional to A'(x) must vanish. This places a condition on the energy function M{E): 



= tr.(AA(i.))./f- 



Af{E) 



(2.16) 



Here /3 = l/T is the inverse temperature and is the chemical potential. This consistency condition imposes one 
relation on the parameters describing the twisted kink solution. With this consistency condition imposed, the general 
inhomogeneous condensate A (a;) satisfying the NLSE (|2.15p solves the gap equation (|2.9p . 

Given this exact solution A(a:) to the gap equation (|2.9p . it is also possible to find the exact single-particle solutions 
to the BdG equation (|2.6p . Furthermore, the diagonal resolvent R{x;E) in (|2.14p encodes all the relevant spectral 
information. For example, the local density of states for fermions in the presence of the condensate is given by 

p{x] E) = ilmtrc {R{x] E + it)) (2.17) 

TT 

where the matrix trace of the resolvent follows trivially from the ansatz (|2.14p : 

tYD {R{x; E)) = 2Af{E) {a{E) + \A{x)\^) (2.18) 

Given the density of states p{E) = J dxp{x]E), all relevant thermodynamic quantities, at finite temperature and 
chemical potential, can be derived from the grand canonical potential 

vI/[A(x);r,/i]=-i pJEp{E)\n(l + e->''^^-^^)+,^\j\x\A{x)\'' . (2.19) 

Since we know p{E) exactly, we can analyze the thermodynamical properties of this model precisely. 
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A. Twisted Kink Crystal Condensate 

The general solution to the NLSE (|2.15[) describes a crystalline condensate [isl. [T^. It is a periodic array of kinks 
that also rotate in the chiral plane, as illustrated in Figure [T] The single chirally-twisted kink was originally found by 
Shei pH using inverse scattering techniques, and subsequently studied in a resolvent approach by Fcinberg and Zee 
[l^. The periodic array of such twisted kinks can be expressed in terms of the elliptic functions: 

A(x) = -A e^'^- A ^1^^^ + -KOa W2) ^""^ ^'^^ ^^'^^'^'^ ^ ' '"'^'^^'^^^ ^ ' ^"^^'^^ ^^'^^^ 
where sc=sn/cn, nd=l/dn are Jacobi elliptic functions, and the functions cr and C are the Weierstrass sigma and zeta 



FIG. 1; The twisted kink crystal condensate of (|2.20|l . shown as the red curve. The blue skeleton surface is shown just to 
illustrate the periodic amplitude modulation and phase winding. 

functions [ssj . chosen to have real and imaginary half-periods; cji ~ K(i'), and = iK' = iK(l — v\ Both periods 
are therefore controlled by a single [real] elliptic parameter < < 1. Note that 773 ~ C(*K') is pure imaginary. The 
parameter A sets the overall scale of the condensate, and 1/A sets the length scale of the crystal. Later, wc will use 
units in which the vacuum mass of the fcrmion is 1, so that A sets the scale relative to the vacuum fermion mass. The 
angular parameter d takes values in the range Q € [0, 4K'(i/)]. The [real] constant A is a fimction of 6 and the elliptic 
parameter v. 

A = A{e, v) = -2i sc(i6l/4; v) nd(i6'/4; v) (2.21) 

For brevity we will usually suppress the explicit dependence of the elliptic functions on the elliptic parameter v. The 
final parameter g is a phase parameter that affects the chiral angle through which the condensate rotates over one 
period L = ^ : 

A(x + L) = e2'^A(x) ; ^ = K (^-z ( W2) + * ns(0/2) - ^ + (2.22) 

where r/ = C(K) is real. Thus the general solution is specified by four real parameters: a scale parameter A, a phase 
parameter q, an angular parameter 0, and the elliptic parameter v. These parameters also parametrize the energy 
spectrum of fermions in such a condensate background, which has two gaps, with band edges Ei < E2 < < E4 as 
shown in the first plot of Figure [H 

El ^ q-X 

E2 - g + A(-l + 2nc2(z6l/4)) 
£^3 = g + A(-l + 2nd^(i6i/4)) 

Ei = q + X (2.23) 

Thus, in terms of the single particle fermion spectrum, the role of the four parameters is as follows: A determines the 
overall energy scale; q determines the overall offset; while 6 and v determine the location and width of the band that 
lies in the gap between the "outer" edges Ei and E4. The simple linear dependence of the energy spectrum on the 
parameters A and g is a direct consequence of the form of the Hamiltonian (j2.5p , and reflects the an important scale 
and shift symmetry described in detail in Section [llll 

For the twisted kink crystal solution (|2.20p . the coefficients a{E) and b{E) in the NLSE (|2.15p are simple polynomials 
of E, with coefficients determined by the band edges: 



aiE) = 2E'-\^E,jE+-^\^E,j - ^(i;, - i?,)^ | (2.24) 
b{E) = 2E-lyEA (2.25) 
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q + X 



-E3- 
E7 



E3 



E\ = q-X 



FIG. 2: The form of the single-particle fermion spectra for the general twisted kink crystal [first figure], showing the central 
value E = q, and the band edges Ej, for j = 1 . . . 4. The second figure shows the special case of the spiral condensate, for which 
the bound band merges with one of the continua. The third figure shows the spectrum for another special case, the real kink 
crystal, which has a charge conjugation symmetry, implying that the offset is g = 0, and the spectrum is symmetric about 0. 
The position of the bands within the gap, and their width, are controlled by the parameters 6 and v. 



Furthermore, the energy function J\f{E) appearing in the resolvent ansatz (|2.14|) also has a very simple form in terms 
of the band edges: 



(2.26) 



Thus, we have an explicit exact expression for the density of states of fcrmions in the presence of such a twisted kink 
condensate field, following from the trace of the resolvent. Within the bands: 



1 aiE) + X^Z 



Here we have defined the function Z{0, v) in terms of the normalized average of |A(x)p over one period: 



(2.27) 



(2.28) 



with V being the Weierstrass V function. Thus, the density of states p{E) is an explicitly known function of the 
energy E, depending parametrically on the four parameters X. q. and v that characterize the solution (|2.20p to 
the gap equation. This parametric dependence enters through the band edge energies Ej in (j2.23p . and through the 
function Z defined in 



B. Spiral Condensate 



An important special case of the general solution (|2.20p is the degenerate case when the bound band of the fermion 
spectrum shrinks and merges with the upper or lower continuum, so that the spectrum has just a single gap, as shown 
in the second plot in Figure [2l This occurs when the angular parameter takes values at its extreme limits: 6* = 
[which implies that E2 = E3 = E4, so that the bound band merges with the upper continuum], 01 6 = 4K' [which 
implies that Ei = E2 = E3, so that the bound band merges with the lower continuum]. The general twisted kink 
crystal condensate (|2.20p reduces to a single plane wave 



(2.29) 



which is clearly a solution to the NLSE (|2.15p . For this condensate the amplitude is constant, while the phase rotates 
at a constant rate, set by q, as shown in Figure [3] The fermion energy spectrum has just one gap, of width 2A, 
centered at q; that is, the band edges lie at Ei = q — X, and E4 = q + X. Correspondingly, the resolvent trace has a 
simplified form, and the spectral function within the continuum bands is simply: 



1 



piE) = - 



\E-q\ 



^^X^-{E-qf 



(2.30) 



FIG. 3: The spiral condensate of (|2.29|l . shown as the red curve. The blue skeleton surface is shown just to illustrate the 
periodic phase winding. In contrast to the twisted kink crystal in Figure [1] for the spiral, the amplitude is constant. 



which we recognize as the spectral function of a constant condensate A = A, shifted in energy by q. 

C. Real Kink Crystal 

Another important special case of the general solution (|2.20p is the case where the condensate is real [relevant for 
the GN2 model], which implies that the BdG Hamiltonian H in (|2.5p has a charge-conjugation symmetry, {H, (T2} = 0, 
which in turn implies that the fermionic spectrum is symmetric, as shown in the third plot of Figure [2j The band 
edges reduce to 

El = — A = — i?4 

E, = -x( l-^] = -E, (2.31) 



,1 + V^/ 

The phase parameter q = 0, and further the angular parameter 9 takes its midpoint value 9 ~ 2K'(zy). Thus the real 




FIG. 4: The real kink crystal condensate of (|2.32|l . shown as the red curve. The blue skeleton surface is shown just to illustrate 
the periodic amplitude modulation and phase winding. For this real kink crystal, the amplitude vanishes each period, and the 
kink rotates through tt [i.e., changes sign] each period. 

kink crystal is described by just two parameters, the scale A and the elliptic parameter v. 

, sn(Ax;^)cn(Ax;^) , _ iy^ 

= tt; — ; ^ = 7; 1=^ ■ (2.32) 

dn(Aa;;z^) {1 + 

The second form of A(a;) in ()2.32p is obtained from the first form by a Landcn transformation [ssj . Over one period, 
L = HJiM^ the condensate changes sign [that is, it rotates through an angle 2(p = —tt], as shown in Figure S) This 
change of sign corresponds to the discrete chiral symmetry of the GN2 model, while the phase rotation (|2.22p of the 
general kink crystal condensate (j2.20p is associated with the continuous chiral symmetry of the NJL2 model. The 
real kink crystal describes the inhomogeneous condensate of the crystalline phase of the GN2 model [iB], and its 
thermodynamics will be discussed below in Section IVIl 



III. THE SCALE AND PHASE SYMMETRY IN NJL2 

In this Section we describe a simple but important symmetry property of the Bogoliubov-de Gennes equation 
that has important consequences for the thermodynamical analysis. The Bogoliubov-de Gennes equation (|2.6p admits 
a family of solutions obtained by rescaling and phase shifting (i.e. making a linear local chiral rotation) a given 
solution: 

A(a;) AA(Ax)e2*«^ 
V'(x) ^ e'9^'^^A^/2^(Ax) (3.1) 



7 



which generates all the linear transformations acting on the energy spectrum: 

E-^XE + q . (3.2) 
In terms of the density of states, the effect of the transformation is: 

piE) - P . (3.3) 

The important physical implication of this symmetry is that when minimizing the grand potential (|2.19p with respect 
to the four parameters A, q, 6 and the minimization with respect to A and q can be done first. If the grand 
potential did not require renormalization, then the minimization with respect to A and q would be trivial. In fact, 
we will show in the next section that even taking into account the renormalization, these symmetries greatly simplify 
the minimization with respect to A and q. 

It is useful to define the "unsealed" and "unshifted" spectrum to be the one with A = 1 and q = 0, so that Ei = — 1, 
and E4 = 1 (in units where the vacuum fcrmion mass is 1). All other spectral functions can be generated from this 
basic solution using the simple transformation (j3.3p . The corresponding density of states will be written as 



[2E^ - {E2 + E:i)E - {E3 - E2y/4: -1 + Z 
' {E^ ~ 1){E - E2){E - E3) 



m-^^ > (3-4) 



where Z = Z{0, v) is defined in ^pJM) . and 



E2 = -l + 2nc2(i6i/4;i/) 

E3 = -l + 2nA^{i6/A;v) (3.5) 



Importantly, p{E) depends parametrically only on the two remaining parameters, 9 and v. This separation of para- 
metric dependences has important consequences for the minimization of the thermodynamic grand potential ()2.19p 
with respect to the parameters. 



A. Transformation Properties of Thermodynamic Quantities 

1. The Grand Potential * 

We begin our discussion with the grand canonical potential \I'[A(a;); T, /^] for the unscaled/unshifted condensate 
A(a;), obtained from ()2.20p by setting the scale parameter A = 1, and the phase parameter q = {]. The grand potential 
is formally divergent in the UV region and has to be renormalizcd, as is well known [lol . [l6| . At finite density and 
nonzero temperature, it is convenient to separate the single particle contribution as 

--/ dEp{E)\n{l + e-^^^-^'^)= dE p{E){E - p.) ~ - dE p{E)\n{l + er'^^^-t'^) (3.6) 

P J-00 JErai^ P J~00 

where i?min = —A/2 — Z/A + . . ., in terms of the momentum cutoff A/2. Only the first term, the zero temperature 
expression, in (j3.6p is divergent. We isolate the divergent terms using the large E behaviour of the density of states 

p{E) « 1 + + . . . (3.7) 



Z 

2E 

The divergent part is 



^ A? Kp Z , , , , 

*div = -^-^-:^lnA . 3.8 
Stt 2tt 27r 

The quadratically and linearly divergent terms are absorbed by definition of the renormalized energy and baryon 
number densities, and the logarithmically divergent term is canceled by the double counting correction 

^ ^ ^ |A(a;)pdx = ^InA (3.9) 



2Nfg^LJo ' ' " 27r 
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where we have used vacuum gap equation jfj^ = InA. Hence the finite renormaUzed grand canonical potential is 

poo 



*,e„[A(.T);r,M] 



dE p(E){E - i^L) - 



dEp{E) ln(l + e-Z^I^-A^I) + ^ + ^ + ^lnA 

OTT ZTT /TT 



(3.10) 



Now we can analyse the effect of the transformation p.3p on the renormaUzed grand canonical potential for the general 
condensate 



A(x) = AA(Ax) e^*"^ 



(3.11) 



The finite temperature (f.t.) contribution [the 2nd term on the r.h.s. of Eq. (|3.10|) ] has the following simple scaling 
behaviour, 



f.t. 



1 

A2 



dEp 



E-q 



^ln(l + e-^l^-^l) 



J- 



dE p{E) ln(l + e 



~l3\E-ii\ 



with the rcscaled variables 



/3 



A2 *,,„[A(x);f,A] 



p.-q 



f.t. 



(3.12) 



^ = A/3 
T 



For the zero temperature contribution (z.t.) in (j3.10p . we start from the expression, 



*rcn[AA(Ax)e2"?^T,^ 



dEp 



E-q 



{E-p + — + -L + 

Stt zn ZTT 



InA 



(3.13) 



(3.14) 



where E^^^^^ = —A/2 — X^Z/A. Here, due to the regularization, the scaling relation analoguous to Eq. p.l2p develops 
anomalous terms akin to the chiral U(l) and scale anomalies, 

4're„[AA(Ax)e2^«^T,/.] = A^ vl/^,,[A(x); T, A] 



H A^ InA + A^^- — 

z.t. 27r 27r 27r 



(3.15) 



Being an UV effect, the extra terms are independent of temperature. Combining Eqs. (|3.12p and (|3.15p . we see that 
the renormalized grand potential for the general condensate in (j3.1ip . is 



vI/„„[A A(Aa;) e^''^-;T,p] = [ M-^en + :flnA+^)-^ 



with the shorthand notation 



2tt 



*„„[A(x);f,A] 



27r 



277 



(3.16) 



(3.17) 



For the sake of compactness in the notation, we will drop the subscript "rcn" from now on, and work exclusively with 
the physical renormalized thermodynamic quantities. 

The grand canonical potential is related to the density p, [not to be confused with the density of states p{E)l], the 
entropy s, and the free energy u: 



— u — pp — Ts 



(3.18) 



Thus we can obtain expressions for the effect of the scaling and phase shifting transformation on the renormalized p, 
s and u as follows: 
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2. Number Density 

From the basic relation p = — ^ , we write ^ = j ^ j ^^id ^ict on p.l6p to obtain: 

2,19^ 1 //^ ^ ^ 



A^(-T^-T-) + - 
\ O jl A TT TT 

A/5+- (3.19) 

TT 



3. Entropy 

From the basic relation s = — ^ , we write -§f = \ -§f^ a-nd act on p.l6p to obtain: 

= As (3.20) 

4- Free Energy 

The transformation property of the free energy now follows directly from the relation (j3.18p : 

u = ^' + /ip + Ts 

^ \^ (u+^\n\\+\qp+^ (3.21) 

\ ZtT J ZTT 

B. Implications for Minimization of the Grand Potential with respect to the phase parameter q 



The minimization of 5* with respect to q can be transformed into minimization with respect to the chemical 

JL = -LJL 

dq A dfi ^ 



potential, due to the symmetry (|3.3p . We write = — and differentiate ^ in (|3.16p with respect to /t 



d-^ ,2/ 5* AM_x2/^. A 



So the g minimization implies 

TTyO = A (3.23) 

Recalling p.l9p and (|3.13p . this means that after minimizing with respect to the phase parameter g, the (period 
averaged) number density is simply proportional to the chemical potential: 

p = ^ (3.24) 

TT 

This remarkable fact is independent of the form of the (complex) condensate, and simply follows from the transforma- 
tion property (|3.ip of the BdG Hamiltonian and its effect on the renormalizcd grand potential, as reflected in (|3.16p . 
Note, of course, that such a relation between p and /i does not arise in the GN2 model, where the condensate is real 
and there is no phase invariancc parameter q. 



C. Implications for Minimization of the Grand Potential 'I' with respect to the scale parameter A 

From (|3.16p . it follows that ^' depends on the scale A explicitly, and also implicitly though the dependence of 
^1 = \&[A;r, A] on T = T/A, and on A = (m - l)/^- Thus we can write 

— = 2A * + — nA H hAM — + AM-- — 

dX \ 2tt J 2tt \ XJ dT \ X J dp 



10 



= 2A(* + |lnA + | + if.s + iA^) (3.25) 
Since "i/ = u — Ts ~ fip, we can express the minimization condition ^ = in terms of the free energy as: 

^ = -|-|hrA + lAP+lf^ (3.26) 

If we impose also the condition ()3.23p arising from the minimization with respect to the phase parameter q, we obtain 
the condition 

Alternatively, we can express these conditions in terms of the thermodynamic quantities for the general condensate 
A(a;) in (|3.1ip . Without using the condition (j3.24p arising from the q minimization, the A minimization condition 
p.26p can be written as 

---^ + h^P+lT'^ + l{p-'i) (3-28) 
After imposing the condition (|3.24p arising from the q minimization, the last term vanishes and wc obtain 

ZA2 ^1^1 , , 

These conditions must hold for any form of the condensate A(x), and will prove very useful in studying the phase 
diagram of both the NJL2 and GN2 models. 

D. Transformation Property of the Consistency Condition 

The final technical ingredient before studying the thermodynamics is the effect of the transformation (j3.ip on the 
consistency condition (|2.16[) . Note that the consistency condition (|2.16p must be satisfied also at finite T and /i, for 
the gap equation to hold. Thus, the energy trace involves the thermodynamical Fermi factor, as in (j2.16p . As with 
the grand potential, density, entropy and free energy, it is useful to express the consistency condition in terms of the 
condensate A(a;) obtained by setting the scale A = 1, and phase q = 0. All we need to know is the effect of the 
transformation p.2p on M{E) . From the form of (|2.26p it is clear that 

M{E)^^n(^^^ (3.30) 

where 

M{E) ^ (3.31) 
{E^ - 1){E - E2)(E - Ez) 



Hence we can write the consistency condition as 

dE M{E) 



27r 1 + e'3(^-'') 



(3.32) 



Note that this integral is finite, even at T = 0, and no renormalization is required. The effect of this condition is to 
express one of the four parameters A, g, 9 and in terms of the others, in a manner depending on T and /i. 

IV. THERMODYNAMICS OF THE SPIRAL CONDENSATE 

Before studying the general twisted kink crystal condensate, we investigate the thermodynamics of the special case 
of the spiral condensate: 

A(a;) = Ae2'«^ (4.1) 
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For this condensate, A(a;) = 1, (i.e., the vacuum fermion mass in our units), and so the thermodynamics is simply 
that of a constant condensate of unit magnitude. The fermion spectrum is now symmetric about 0, and so we can 
immediately write an expression for the corresponding grand potential ^: 

^ = -— - ^ ^P. , ^ ln((l + e-'^(^-«)(l + e-'3(^^+A))) (4.2) 

The full grand potential ^E* for the spiral condensate (|4.ip is then obtained using p.l6p . Next wc minimize the full 
grand potential with respect to q and A. 



be written explicitly as 



A. Minimization with respect to the phase parameter q 

At T = 0, wc see from (|4.2p that 'I' = — independent of //, so that p = 0. Therefore, the condition (|3.23p . arising 
from the minimization with respect to q, implies that /i = at T = 0. In other words, q — ii, so that the chemical 
potential lies at the center of the gap in the single-particle fermionic spectrum. With q = /z, the spiral condensate 
(|4.1[) is the "chiral spiral" solution proposed in [1^. At nonzero temperature, the q minimization condition (|3.23p can 

fi = np = 77— 

O/J, 

= 2sinh(/3/i) / dE , -. (4.3) 

Jl VE^ - 1 (1 + e'3(-^-A))(l + e/5(^+A)) 

At low temperatures, T ^ 1, the main contribution to the energy integrals in (|4.3p comes from near the upper band 
edge E = 1. So we approximate the density of states as: 

and (14.31) becomes 



A « %/2sinh(/3A)^ di?-^^=e^^e-'^(^-'^)e-'^(^+^) = ^^e-^sinh(/3A) (4.5) 

This also requires /i=0, leading again to the chiral spiral solution with q ~ fi. Indeed, it is easy to verify numerically 
that the finite temperature equation (j4.3p has a solution only at /i = 0, for all temperature T. Another argument in 
favor of fi = q at all temperatures is that instead of minimizing with respect to q, we can minimize with respect to 
/t. Since 'I' is symmetric under /i — > — /i, there must be a stationary point at /i = 0, i.e., fJ- = q. That it is a minimum 
can easily be seen by looking at the sign of the 2nd derivative (Taylor expansion of the integrand). Other minima 
(which could only come in pairs) are ruled out numerically. 

Thus, we conclude that the minimization of the grand potential with respect to the phase parameter q leads to 
q = ^ for all temperature T, so /i always lies at the center of the gap. As should be clear from this discussion, this 
fact can be traced directly to the phase transformation symmetry in (|3.ip . 

Another immediate consequence of g = is that the grand potential for the chiral spiral has a simple /i dependence. 
This follows because p.23p with jl ~Q implies that ^ is independent of the chemical potential /i. Indeed, when fi ^ 0, 
the grand potential (|4.2p can then be written as 

^ _± _ ^ r dE , ^ ln(l + e-^^) (4.6) 

Then the general relation (|3.16p implies that for the chiral spiral condensate the full grand potential is 

^ = ^ (In a2 - 1) - — dE , ^ ln(l + e-^^) - ^ (4.7) 

Thus, the grand potential for the chiral spiral has a simple p, dependence, and it is clear that p = —d'^ /dp = p/ir. 
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B. Minimization with respect to the scale parameter A 

From (|4.7p it also follows that the scale parameter A is determined only by T, independent of the chemical potential 
/i. Indeed, minimizing (j4.7p with respect to A, we obtain the equation for the thermal mass scale A(r): 

= X—-—r dE^= ln(l + e-f'^) + A- H dE , ^ (4.8) 

It is a simple exercise to show that this is equivalent to the general A minimization condition (|3.29p . expressed in 
terms of the entropy and the free energy. Note that this equation does not involve the chemical potential fi, so the 
thermal mass scale A(r) must be independent of /i. 

At r = 0, (|4.8p reduces to A In A = 0, which implies that A(r = 0) = !, and the grand potential is simply 

^T=o^__L_/^ (4 9) 

47r 2tt ^ ' 

For small but nonzero temperature, the scale parameter A receives an exponentially small finite T correction, found 
by approximating the energy integrals in (|4.8p : 

A(T) - 1 - %/2^e"^/^ , T < 1 (4.10) 



Applying the same approximation to the minimized grand potential in ()4.7p we find the leading small T correction to 
the grand potential: 



1 ./^r'^-i/T ,4 11) 



Att 2tt V tt 

For general T, the temperature dependent mass scale A(T) can be obtained numerically from (|4.8[) . The scale A(r) 
decreases monotonically from the value A = 1 at T = 0, and vanishes at a critical temperature 

= — w 0.566933 (4.12) 

TT 

At this temperature, T ~ T^, the system undergoes a phase transition to a massless phase. Interestingly, this phase 
transition is independent of the chemical potential /i, as follows from the fact that A(r) is independent of ^. We 
can trace this fact directly to the simple form (|4.7p of the grand potential for the chiral spiral condensate, after 
minimization with respect to the scale parameter q. 

Just below Tc the dependence of A on T is nonanalytic, as can be seen from the following argument. After integrating 
by parts the second integral in (|4.8p . and expanding the Fermi factor we obtain 

= InA + 2^(-l)"+i / dEj==e~"P^ = InA + 2^(-l)"+i/io("/3) (4.13) 

n— 1 ^ ^ n—l 

where ^^0(2^) is the modified Bessel function. To obtain the critical exponents near the phase transition, wc analyze 
this equation near small values of A. Since /3 = A/T we expand the Bcsscl functions around zero: 



= lnA + 2^(-l)"+i 



-lM^)-7-^ In -1+7-1 



(4.14) 



n=l 

Here 7 is Euler's constant. The n sums can be evaluated in terms of the Riemann zcta function, leading to 



= ln(T.) - 7 - (4.15) 

In particular, at the phase transition where A = 0, the critical temperature is found to be = e'''/7r, and to the 
leading order in (T^ - T), for T < T^: 



/ 2T 

HT) = J _7^,(l2) VTc-T+...^ 3.06 VTciTc - T) + . . . (4.16) 
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Thus, for the spiral condensate (|4.ip . the thermodynamic phase diagram is given in Figure [5l showing the phase 
transition at T = T^,, independent of /x. After minimizing the grand potential we learn that in the region T < T^, 
the pitch angle q of the spiral condensate is directly proportional to the chemical potential, q = /i, independent of T, 
while the amplitude A(r) is just a function of temperature [vanishing at Tc], independent of /i. 




FIG. 5: The phase diagram of the NJL2 model. The tricritical point is marked at /itc ~ and Ttc ~ e^/yr « 0.5669. Below Tc 
the condensate has the form of the spiral condensate (|4.ip . with q = n- 



V. THERMODYNAMICS OF THE TWISTED CRYSTAL CONDENSATE 



In this section we develop results for the thermodynamics of the NJL2 model with the twisted kink crystal conden- 
sate, that is the general solution of the inhomogeneous gap equation. Recall that the twisted kink crystal condensate 
(|2.20p is characterized by 4 parameters: the scale parameter A, the phase parameter 5, an angular parameter 6, and 
the elliptic parameter v. We first consider the situation analytically at T = 0, then at nonzero T. 



A. Twisted Kink Crystal at T = 

At T = there are some significant simplifications. First of all, the Fermi factor becomes a step function that acts 
as a cutoff of the energy integrals. Thus, in the consistency condition (|3.32p . we use 

Q{fi-E) (5.1) 



1 + e'^(^-^) 



If we assume that /t is in the upper gap, then both the lower continuum and the bound band are completely filled. 
Thus the consistency condition p.32p reads (for details, see [l4| ) 



-1 rE3 ^ 



- / dEN{E) + / dEN{E) = K' (5.2) 



Therefore, the consistency condition forces 6 = 4K', which is precisely the spiral condensate case. In this limit, the 
band shrinks and joins the negative energy continuum, leaving just the single-gap spectrum of the spiral condensate. 
Then the q minimization leads to g = as described in the previous section, and we find the preferred condensate 
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to be the chiral spiral. Similarly, if we assume that /t is in the lower gap, then we get = 0, that is also the spiral 
condensate limit; in this limit the bound band joins to the positive energy continuum. Once again, minimization with 
respect to the phase parameter q leads to g = /i, so the condensate is the chiral spiral. 

The only other possibility is fl lying inside the bound band. In this case the consistency condition (|2.16p leads to 
an expression for the Fermi energy: 

Ef ^ Ef{0,v) =TLc{ie/2\v) (5.3) 



l.U 




FIG. 6: Plot of TTp{^,u) [light blue surface] and _Ef(C,J^) [darker blue surface], as functions of u and ^, where 6 = 4K'^. The 
surfaces intersect at ^ = 1/2, which means 6 = 2K'. 



On the other hand, at T = 0, this Fermi energy is simply the chemical potential. Minimization with respect to 
the phase parameter q leads to the relation fi = it p. We can evaluate the density obtained by filling up to the Fermi 
energy Ef{0,p): 



1 



1 



cii{ie/2]v) - dn{ie/2;v) 
{ie/2-v) \cu{ie /2]v) + An{ie /2]v) 



(5.4) 



The simultaneous solution of these two conditions, namely Ef = Trp, has the unique solution 9 = 2K', for all ly, as can 
be seen from Figure [S) Evaluating the free energy for this solution we obtain the following function of the remaining 
parameters A and v. 



2A2 



^(l + v^)2 



EM 



In 



2A 



1 + 



27r 



(5.5) 



Minimizing £(A, v) with respect to v, we find that we are forced to = 1, which means 



(5.6) 



from which we recognize the T = grand potential (|4.7p of the chiral spiral solution. Thus, once again, the mini- 
mization forces us to the chiral spiral condensate solution, at T = 0. 



B. Twisted Kink Crystal at < T < 1 



The minimization of the grand potential at T = shows that the preferred twisted kink crystal configuration is 
the chiral spiral, with the chemical potential sitting in the middle of the gap (i.e = [p, — q)/\ = 0). For this 
solution, the angular parameter 6 takes the values or 4K'. We will consider the latter case [a similar argument 
applies for the other choice]. Now consider the stability of this chiral spiral for T nonzero but small. If we change 
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T slightly away from 0, then the consistency condition that sets the angular parameter 6 = 4K', will also change 
slightly, and the preferred value of 9 will shift away from 4K'. We write 9 = 4K' — 4e, where e ^ T <C 1. This 
changes the single-particle spectrum by producing a very narrow band very close to the lower band edge Ei ~ — 1. 
With 9 = 4K' — 4e, the band edges p.5p and the averaged amplitude (|2.28p take the form: 

M^,'^) ~ -l + 2jye^ + ... 
hid, I') « -l + 2e^ + ... 

Z{9,v) « l-h{v)€^ + ... (5.7) 

where h{v) = 2v — 2 + 4E(i/)/K(i^). We now calculate the small T correction to the grand potential \E' for this 
configuration. As usual, we split the grand potential into zero temperature and finite temperature parts. For small 
T, we use the fact that ln(l + g-l'^^^"'^)!) « e'l'^^-^-'^)! « e"'^!^! (recall that /i = for T = 0). 

/ — I pEz I p — 1 . poo _ pEs . \ 

dEp{E){E- / dEp{E){E - ft) - T [ / dE p{E)e'^^ + / dE p{E)e'^'^ + / d£;p(£:)e^^ 

-OO JB2 \J -00 Jl Je2 J 

) poo r^^^ \ 

= < ^ >T=o -Te-l^ { I dx p{-x - l)e'^'' + dx p{x + l)e-^'' + dx p{x - l)e^^'' ] (5.8) 



In the small T limit, the continuum integrals are dominated by the region x ^ (i.e near the band edges). The 
spectral function around the band edges has the behavior: 

Ft' 1 

00 / rp 

dx p{x + l)e-^'' « \ — + 0{e^) (5.10) 

V ZTT 

Note that the lower contimmm leads to an 0(e) correction, while the upper continuum leads to an O(e^) correction. 
The band integral also leads to an 0(e) correction. This can be seen by changing the integration variable x e^x: 

:'dEpix-l)e^^ ^efdx e^-' -7(1 + -)- M^)/^ ^ , ^ . - 2E(.)/K(.) ^ ^ ^^^^ ^^^^^^ 

2i/£2 J21/ IT ^ x{x — 2v){2 — x) J21/ TT ^J x{x — 2v){2 — x) 

An important observation is that this function L{v) < for all v G [0, 1]. 

Finally, we use the general transformation of the grand potential (|3.16p to deduce the grand potential for the twisted 
kink crystal. Since after minimization, ipT=o = — l/(47r), the small T correction to A does not contribute to the full 
grand potential \l/ (as in (j4.1ip ). and the full grand potential is found to be: 



The first term in parentheses is just the low T grand potential for the chiral spiral, as in (j4.1ip . Since Liv) is always 
negative, and K(y^) is always positive, we see that the system is unstable with respect to the opening of a gap near 
the lower continuum edge. In other words, at small T the minimization of the grand potential reduces the general 
twisted kink crystal condensate to the chiral spiral condensate, just as at T = 0. 



C. Numerical results for the thermodynamics of the twisted kink crystal condensate 

The previous two sections have shown that at T = and for small T, the chiral spiral condensate is the thcrmody- 
namically preferred form of the more general twisted kink crystal condensate. We have also checked this conclusion 
numerically at various locations on the phase diagram, and we find that throughout the phase diagram the chiral 
spiral is the thermodynamically preferred limit of the general twisted kink crystal solution of the inhomogeneous gap 
equation. 
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VI. THERMODYNAMICS OF THE REAL KINK CRYSTAL AND THE GNa MODEL 



The phase diagram of the GN2 model, which has just a discrete chiral symmetry instead of the continuous chiral 
symmetry of the NJL2 model, is by now well understood [T6l.[3^. But we revisit it here briefly, with a new perspective. 
The analysis of this paper shows that the key to understanding the phase diagram of the NJL2 model is the behavior 
(|3.16p of the renormalized grand potential under the rescaling and shifting transformation (|3.ip . But in the GN2 
model there is no pseudoscalar interaction, so the condensate is real. Thus, there is no symmetry corresponding to 
a phase rotation of the condensate. In other words, q = 0. The angular parameter 9 in the solution (|2.20p of the 
inhomogeneous gap equation is also zero, as the condensate cannot wind by an arbitrary phase as it goes through one 
period. Furthermore, there is no need to impose any consistency condition on the solution of the gap equation: since 
the pseudoscalar condensate 11 is identically zero, there is no condition arising from its variation, which means that 
the off-diagonal terms in (|2.1ip play no role. Thus, the general solution (|2.20p of the inhomogeneous gap equation 



o.s-i 




FIG. 7: Phase diagram of the GN2 model. The tricritical point is at /itc = 0.608 and Ttc = 0.318. In the region of /i > 2/-K, 
the massless and massive phases are separated by a crystalline phase. 

simplifies to the real kink crystal solution in (|2.32p . which depends on just two parameters, the scale A and the elliptic 
parameter v. Concerning the grand potential, the key formula is now p.l6p . with q set to 0: 

*„„[AA(A.T);r,/i] = A2 (^*,en[A(a;);T/A,A*/A] + ^lnA^ (6.1) 

The last term reflects the anomalous behavior of the grand potential under the rescaling of the condensate by A. 



A. Real kink crystal at T = 

From previous work [l^ , we know explicit expressions for the thermodynamical quantities at T = as functions of 
the elliptic parameter u. For A = 1, the density is 

1 



The free energy is 



(6.2) 



Z{D) 1 E(^) 

— — (In - 1) + — — — — (6.3) 
47r 27r K(i/) 
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where 

which is just Z{Q^v) in (|2.28p . evaluated at = 1K!(y). The function E(i>) is the complete elliptic integral of the 
second kind [s^ . Thus we can write the T = grand potential as 

* = £ - iip 

= A2(/i(i>) + /2(i>) lnA)-/.A/3(i>) (6.5) 



where 



ZTT 

^^^^ - 21^ 
Minimizing ^ with respect to A and v leads to two equations: 

^ = 2A (^/i + i /2 + /2 InA^ - /i/3 = 

^ = A{A(/{+/^lnA)-M/a = (6.9) 

Simultaneous solution of these conditions leads to a complicated- looking expression for In A, that actually simplifies 
dramatically; 



In A 



/3/{-2/^(/i + /2/2) 
2/2 /a - /s /2 

= -iln:> (6.10) 



In showing this remarkable reduction we use the property 

dZ {Z-vf 



dv 4i>(l - v) 

Inserting this result for A back into the minimization conditions (|6.9p wc find the minimized values at T = 0: 



m = 4^ (6.13) 

The critical value of chemical potential, /ic = f , corresponding to the baryon mass [1, [i3j[i3> obtained at i> = 1, 
in agreement with known results p^ . 

B. Real kink crystal at T < 1 

At nonzero temperature, the minimization with respect to A and v leads to T dependent expressions for the chemical 
potential and the scale factor A, as functions of the elliptic parameter i>, generalizing the T = expressions (|6.121 
I6.13p . This can be done numerically, as in (l6j . but here we find analytic expressions valid in the small T limit. First, 
we note that the T = chemical potential in (|6.12[) lies in the upper gap [sec Figure [8]: 



E-i(y) = A(i>)\/T~~^ < m(i>) < E^{v) ^ A(£>) (6.14) 
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FIG. 8: Plot of the chemical potential [center line], and the band edge energies, as a function of the elliptic parameter v. Note 
that /i(^) lies in the gap for all v, and moreover, it is slightly closer to the upper band edge, £4, than to the lower band edge, 
E3. 



Since /i is in the gap, at small T there is an exponentially small factor in the corrections to thermodynamic quantities 
going like 



exp[— l/i — nearest band edge|/r] 
Furthermore, for all fi is closer to E4 than to E^. Thus, we can write as a leading approximation 



•9 = -T j dEp{E) In [l + e 

JEa 



0{E-Ei) 



We expand the spectral function in the vicinity of the nearer band edge, E4: 

2E^- - {Ei + El) + Z 



P{E) 



2^^{E^-El){E^-El) 
1 Ej-Ei + X^Z 1 



27r y/2Ei{El - El) VE - E^ 



Oi^E - Ei 



Thus, 



where 



1 



E(^) 
K(^) 



(6.15) 

(6.16) 
(6.17) 
(6.18) 

(6.19) 
(6.20) 

(6.21) 
(6.22) 



We now minimize 5* with respect to A and keeping the leading small T corrections to the T = results of the 
previous section. We find, after some straightforward algebra. 



p{y,T) 



2E(i>) /2T (1 - i))'K{v) 



cxp 



2E(£>) 



X/ N 1 h^T 1 
A(.,T) ^ ^^__exp 



2E(z>) 



(6.23) 
(6.24) 
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These small T corrections are plotted in Figure [9l and are in very good agreement with the numerical results found 
in [l^ , and plotted in Figure [71 Already these corrections indicate the existence of a crystalline phase in which 
the condensate scale A and the period (set by the elliptic parameter i/) are dependent on both T and fi. This is 
in contrast to the phase diagram of the NJL2 model, shown in Figure [5l where the phase transition line is only a 
function of T, and the scale parameter A is independent of the chemical potential With this perspective we can 
trace this fundamental difference in the phase diagrams directly to the fundamental difference between the discrete 
and continous chiral symmetry of the two models. 



FIG. 9: Plots of the chemical potential, as a function of T, for various values of the elliptic parameter D, keeping just the 
leading small T behavior, as in (|04ll . This plot is in remarkable agreement with the numerical results shown in Figure [7] 



We complete our analysis of the phase diagram of the NJL2 and GN2 models by analyzing another region of 
the phase diagram, using the Ginzburg-Landau expansion of the grand potential ^P. In the previous sections we 
obtained analytic results at and near T = 0, but the Ginzburg-Landau approach permits a certain degree of analytic 
information about another region of the phase diagram, in the vicinity of the tricritical point. Expanding in powers 
of the condensate and its derivatives, the renormalized grand potential density may be expressed as: 



0.5 




T 



2 



VII. GINZBURG-LANDAU ANALYSIS 



*GL = ao + a2| A|2 + asim (AA'*) + a4(| A]"* 



+ |Af ) + a5lm((A"-3|A|2 



A)A'*) 



+a6(2|A|<' + SlAplA'p + 2ReA'2A*2 + |A'f ) + . . . 
The coefficients a„(r, /z) are the following functions of T and n [l6l |: 



(7.1) 
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1 



2127^537-4 



Reip 



(4) 



1 .Pf, 

2+^2^ 



(7.2) 



Keeping terms up to a certain order in this expansion, and inserting into the gap equation (|2.9p . we obtain an 
equation (the Ginzburg-Landau equation) for the condensate A. Remarkably, for the NJL2 and GN2 models, this 
hierarchy of equations can be solved to all orders H^l- If we expand up to a2, then the Ginzburg-Landau (GL) 

equation is simply A = 0, so we learn nothing about the phase diagram. To this order the system appears to be just 
a free masslcss Fermi gas. If we expand ^' up to 03, then the Ginzburg-Landau (GL) equation reads: 



A' 



."2 . 

I — A 

as 







A = A cxp 



.a2 

I — X 

as 



(7.3) 



This has the form of the spiral condensate studied in Section IIIBI This spiral condensate has constant magnitude, 
|Ap = A^, and also ^ (AA*' — A* A') = so that when we evaluate the grand potential on this solution, we 

find 



GL 



"0 



(7.4) 



So the grand potential is independent of A, and this is again no different from a free massless phase. If we expand up 
to Q!4, then the Ginzburg-Landau (GL) equation is the NLSE equation: 



(A" - 2|A|2A) - z^A' - ^A = 



(7.5) 



0:4 

The general bounded solution of this equation is the twisted kink crystal described in Section [II Al 

The general pattern is the following: to order a^, the GL equation is a differential equation of order (fc — 2), and 
the general solution corresponds to a finite gap Dirac problem with {k — 2) gaps, or (k — 1) bands (including the 
positive and negative continuum bands). For example, the a2 equation led to A = 0, which is the free system with no 

gaps. The as equation leads to A = A exp x , which has precisely one gap. The a^ equation, the NLSE (|2.15p . 
has as its general solution a system with two gaps, as shown in the first plot of Figure [21 In general the solution with 
{k — 2) gaps requires 2(fc — 2) parameters for the solution, and these parameters can be thought of as labelling the 
band edges. Let us write 



GL 



1=2 



Then the GL equation to order k is: 



E 

1=2 



ai(T,fi) 



SJi[A,A\A" 
5A*{x) 



^ 



(7.6) 



(7.7) 



These equations define the AKNS hierarchy, and (|7.7p is also known as the Novikov equation. Formal expressions 
exist for their solution in terms of multi-dimensional theta functions [26l . [27j , although these are cumbersome to work 
with beyond the twisted kink crystal solution. The remarkable integrability properties of the AKNS hierarchy implies 
that the solution to the NLSE satisfies the Novikov equations to all orders, with suitable choices of parameters, as 
shown also in UM. 



A. Ginzburg-Landau expansion for the GN2 model 



It is instructive to see how successive orders of the Ginzburg-Landau expansion reveal more and more about the 
exact phase diagram. In the GN2 system, the condensate is real, and we write it as A = 0, and all odd-index terms 
of the Ginzburg-Landau expansion vanish. The Dirac spectrum is now symmetric about 0, and so the band edges of 
the finite-gap solutions come in ± pairs, as in the 2-gap case depicted in the third plot of Figure [2j Therefore, we 
only need half as many parameters at a given order to describe the solution. The grand potential density simplifies to 



'I'gl = ao + a2(f>'^ + a4 [ (p'^ + (f)''^ - - 



24,^ 1002 



- 



(7.8) 



The tricritical point is defined as the point where the first two nontrivial coefficients, a2{T,fi) and a4{T,fi) vanish: 
a2{T, n) = ai{T, /i) = ^ = 0.318329 , ^tc = 0.608221 (7.9) 
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1, Ginzburg-Landau expansion to 0(ai) for the GN2 model 



As mentioned above, the expansion to 0{a2) yields no information. The next nontrivial order, to 0(04), leads to 
the following GL equation, which is a special case of the NLSE (|7.5p 



2(f)* 10 = 

Q4 



The general solution can be written as 

with the identification of the scale parameter A as 

A2 = 



+ (1 + = 



a2 



(7.10) 



(7.11) 



(7.12) 



Notice that in the GL approach, we get explicit expressions for the dependence of the solution's parameters in terms 
of T and /i. An important comment is that since > 0, this expression tells us that this inhomogcneous solution 



only makes sense in regions of the (T, /i) plane where 
solution in (mT|) 



Q2 
"4 



> 0. Using the following identities satisfied by the 



(02)" ^ 60'^-4(l + z/)A2 02 + 2i.A4 
we can write the grand potential density to this order as 

^GL = 00+020^+^4 Q(l + i^)A2 02 + ii^A''^ 



OO + g02 ' 



(7.13) 
(7.14) 
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Here we have used the above expression (|7.12l) for A^. Averaging over one period, we use 
and again using (|7.12p we find 



(7.15) 

A2(l-E(^)/K(^)), 
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(7.16) 



Note that the function F{v) is a smooth function interpolating monotonically between F{Q) ~ and i^(l) = 1. We 
have written (^)ql''*^' like this in order to compare with the homogeneous ansatz: = A. Then 

^ homogeneous 



GL 



oq + a2A + 04A 



Minimizing with respect to A^, we obtain the condition A^ = —02/(204), and at this minimum 



/.T.\ homogeneous 
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Therefore, we can write 



= Oo 



/^\ homogeneous /^\ crystal 

\ / GL \ / GL 
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(7.17) 



(7.18) 



(7.19) 



An important observation is that at the values oi v = 1 and v = Q, the minimized grand potential reduces to that of 
the homogeneous and the massless condensates (recall E/K(i^ = 1) = 0,E/K(i^ = 0) = 1): 



(*)-r^v = 1) = «o 
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GL 



(7.20) 
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FIG. 10: The phase diagram of the GN2 model, based on a Ginzburg-Landau expansion to the lowest nontrivial order: 0(a4). 
The blue region is the region in which the massive homogeneous condensate (j!> = A has a lower grand potential. The white region 
is the region in which the homogeneous massless condensate = has a lower grand potential, or where only the massless 
condensate exists, because in (|7.12p is negative. These regions meet at the tricritical point: Ttc ~ 0.318, /itc ~ 0.608. 

This behavior is depicted in Figure [TOl where the grand potential of the crystal condensate lies between that of the 
massless and massive homogeneous phases, interpolating between them as a function of v. Minimizing with respect 
to V pushes us to the massive homogeneous phase in the blue region, but to the massless homogeneous phase in the 
white region. Thus, at this order of the GL expansion, even though the solution to the GL equation has the form of 
a crystalline condensate, the thermodynamic minimum is a constant condensate, cither zero or non-zero, but always 
constant. We show in the next section that this picture changes significantly at the next order. 



2. Ginzburg-Landau expansion to Ofaa) for the GN2 model 

Going to the next non-trivial order beyond the level defining the tricritical point, we expand the grand potential 
density in powers of the real condensate field cj) and its derivatives [we drop the total derivative terms as these are not 
important for this argument]: 



The GL equation is now a fourth-order equation: 



U"" - 1002^" - 100(0')' + G^') + — i-r + 203) + 11^0 = 

The simplest solution is a homogeneous condensate, = A, with massless and massive solutions 

A = (massless homogeneous phase) 
a2 



0L2 



— = ^ xl^-^ 



1± Wl- 



(7.21) 



(7.22) 



(massive homogeneous phase) (7.23) 



3^6 ^u:g ^<-*ti \ V '-^4 

The general solution to (|7.22p is very complicated, but we can use the inhomogeneous solution to the NLSE 

=: A^/^sn(Aa;,^^) (7.24) 
A similar idea was used in an analogous condensed matter model in [s^ . This solution satisfies the nonlinear equations: 

-0" + 203 = (l + jy)A20 



100V' - 100(0')' +60^) = (v^ + 



(7.25) 
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FIG. 11: The phase diagram of the GN2 model, based on a Ginzburg-Landau expansion to the lowest nontrivial order: ©(ae)- 
The blue region is the region in which the massive homogeneous condensate (f) = X has the lowest grand potential. The red 
region is the region in which the crystalline condensate (|7.24p has the lowest grand potential. In the remaining [white] region, 
the homogeneous massless condensate <j) = has the lowest grand potential. Note that the crystal phase region begins at the 
tricritical point: Ttc = 0.318, fitc = 0.608. A close-up of this region is shown in Figure [12] 



Thus, comparing with the GL equation (|7.22p . we see that (f> satisfies the GL equation (|7.22p provided we identify: 

+ + 1) Q,g + + 1) Q,g 

This condition leads to two solutions 



+ . /: + , , .. ^ = (7.26) 



^±-~2(^2+4^ + l)-|l±(^l- (^ + 1)2 ^) 1 (7-27) 



Evaluated on the crystalline solution, the grand potential is 



\^„. I 1 I I ^(11 o.. /I I ,,^ I I "1 o,,2 I c,, I 1 ,.2. I A- I i\ 



(*>GL = ao + A"a2 (^1 - K j ^ ~Y' (^^ + " + j + \ + + I - [v' + + 1)- j (7.28) 

This is just a function of T and /x (through the a's) and the elliptic parameter v, because A is given by the solutions 
in (j7.27|) . We can therefore evaluate the grand potential throughout the (T, /i) plane and ask where it is lower than 
the grand potential of the homogeneous phase. The result is shown in Figure [11] which shows the existence of a 
crystalline phase in a small region in the vicinity of the tricritical point. This is a region in which the grand potential 
of the crystalline condensate is lower than that of the massless or massive homogeneous condensate. On the upper 
edge, V = Q and the scale of the crystalline condensate vanishes as it reduces to a massless phase; on the lower edge, 
J/ = 1, and the period of the crystalline condensate diverges as it reduces to a homogeneous massive phase. The form 
of this region matches very well with the full crystalline region, near the tricritical point, as shown by the close-up 
view in Figure 1121 Going to higher orders of the GL expansion, this crystalline region grows, and eventually covers 
the entire region given by the exact numerics [l6j . 



B. Ginzburg-Landau expansion for the NJL2 model 



In contrast to the GN2 model, the a^, term in (|7.ip is present in the GL expansion of the NJL2 model, as the 
condensate A is complex. The "tricritical" point is defined as the point where the two lowest nontrivial coefficients. 
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FIG. 12: A close-up view of the crystalline region in the phase diagram of the GN2 model, near the tricritical point, based 
on a Ginzburg-Landau expansion to the lowest nontrivial order: 0{ag). The red shaded region is the crystalline region seen 
at this order of the GL expansion, while the solid black lines mark the edges of the true crystalline region found numerically 
from the exact grand potential [1^ . The agreement is excellent near the tricritical point and near the LOFF boundary with 
the massless phase. 



a2 (T, fJ-) and as (T, fi) , vanish: 

a2(r, //) = 03(7, /i) = => Ttc = 0.566 , Htc ^ 



(7.29) 



1. Ginzburg-Landau expansion to 0(a-j,) for the NJL2 model 



The first nontrivial order, to 0{a^), leads to the GL equation 



A'-^^A = 
as 



A = A exp 



.a2 

I — X 



(7.30) 



But for this solution, even though this condensate is crystalline, the grand potential is (vE')gl = ao- Thus, the phase 
diagram is simply that of a massless phase. The only thing we learn at this level of the GL expansion is the existence 
of the tricritical point at T = 0.5669 and /i = 0. This is analogous to the situation of the GL expansion of the 
GN2 model to its first nontrivial order. 0(04), where the solution of the GL equation has a crystalline form, but this 
crystalline condensate does not appear in the phase diagram at that order, as discussed in Section [VII A II 



2. Ginzburg-Landau expansion to O(cti) for the NJL2 model 

Going to the next non-trivial order beyond the level defining the tricritical point, namely to 0(04), wc obtain the 
GL equation of NLSE form in (|7.5p . Adapting the solution in Section III Ap , we can write the general solution as 



which satisfies 



^ - - \(A. + .K0.(.//2) [^^-(-<(^^/2) + ^) + ^^^3/2] 
- A" -f- 2A|Ap = -2igA A' + X^{^W{ie /2) - q^)A 



(7.31) 



(7.32) 
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Identifying the terms with the NLSE equation we deduce q 



a2 ^ 
2a4 4a2Q^4 



and A must satisfy: 
2 



-3V{ie/2) 



(7.33) 



Note that {—V{i9/2)) > 0. Thus, this inhomogeneous crystal condensate only makes sense in regions of the {T,ijl) 

> 0. 



plane where ( — §^ 
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FIG. 13: The phase diagram of the NJL2 model, based on a Ginzburg-Landau expansion to 0{ai). The purple region is the 
region in which the spiral condensate (|4.H) has a lower grand potential. The red region is the region in which the massless 
condensate A = has a lower grand potential, or in which only the massless condensate exists, because in (|7.33|l is negative. 



Now evaluating the averaged potential on this solution, we find 
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(7.34) 



which should be compared with the corresponding GN2 expression (j7.16p . [Indeed, setting 6 ~ 2K', and 03 = 0, we 
recover the GN2 formulas]. We note that < F{iy,9) < 1. We now compare the averaged potential for the crystal 
with that obtained from a spiral ansatz: 



With this spiral ansatz we find 



spiral 



^spiral ^ g22qa: 



oo + A'^q;2 - 2qX'^a3 + (A* + Aq^y)ai 



Minimizing with respect to q we find q = 03/04, and further minimizing with respect to A^ we find 

2o4 4o2Q;4 



(7.35) 
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(7.37) 



which should be compared with (|7.33p . Furthermore, evaluating the averaged potential on this spiral condensate we 
find 
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(7.38) 
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which should be compared with (|7.34p . 

Now we combine this expression with the positivity condition on in (j7.12p to obtain the phase diagram in Figure 
[T3l Just as in the GN2 case, here in the in the NJL2 model, by going one step beyond the first nontrivial order of the GL 
expansion [i.e., one step beyond the order that defines the tricritical point] we see the appearance of a crystalline phase 
in the phase diagram, in the region near the tricritical point. For the NJL2 model, the condensate of this crystalline 
phase, derived from this GL approach to this order, has the form of the chiral spiral after minimization of the 
grand potential. This Ginzburg-Landau analysis confirms once again that the chiral spiral is the thcrmodynamically 
preferred form of the inhomogcneous condensate, in the applicable part of the phase diagram. The pattern is fairly 
clear: going to higher orders of the GL expansion, the crystalline region grows, and eventually covers the entire region 
given by the exact numerics, as shown in Figure O 



VIII. CONCLUSIONS 



We have used the exact crystalline solutions to the inhomogcneous gap equation of the NJL2 model, found in 
[1 .[1^. to probe the thermodynamic phase diagram of the NJL2 and GN2 models at finite density and temperature. 
Using a combination of exact, numerical and Ginzburg-Landau approaches, we have shown that for the NJL2 model 
the thermodynamically preferred condensate in the region T < T^is the helical chiral spiral of [l^ . The same methods 
have been applied to the GN2 model, confirming previous numerical results [l6l |. A key new idea in our analysis is 
the exploitation of the behavior of the grand potential under the rescaling and phase rotation transformations (|3.ip . 
which affect the renormalized grand potential as in (j3.16p . This observation greatly facilitates the minimization of 
the renormalized grand potential with respect to the parameters A and q. We are also able to trace in a very explicit 
manner the consequences for the phase diagram of the fact that the GN2 model has a discrete chiral symmetry, 
while the the NJL2 model has a continuous chiral symmetry. These one dimensional models are somewhat special, 
due to the rich integrability structure underlying their gap equation. So, we studied these models also using the 
Ginzburg-Landau approach, which does not necessarily rely on this integrability structure. We found that in both 
the NJL2 and GN2 models the crystalline region appears at the order of the Ginzburg-Landau expansion one step 
beyond the first nontrivial order, which is used to identify the relevant tricritical point. It would be interesting to 
study this point systematically in higher dimensional models, where the search for crystalline phases is considerably 
more difficult [2Qi . i35. . ,36 . .37j | . It would also be interesting to study the NJL2 system on the lattice, complementing 
the GN2 work of [Ull^l, and the recent Monte Carlo formulations in [3^[3^, ll^]. Other interesting effects include 
the study of an isospin chemical potential [4l| , and going beyond the leading large N approximation [4^ 1 . 
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